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The effort in the ORNL-SNAP shielding program is directed toward the development 
and verification of computer codes using numerical solutions to the transport equation 
for the design of optimized radiation shields for SNAP power systems. A brief discus- 
sion is given for the major areas of the SNAP shielding program, which are cross-section 
development, transport code development, and integral experiments. Detailed results 
are presented for the integral experiments utilizing the TSF-SNAP reactor. Calculated 
results are compared with experiment for neutron and gamma-ray spectra from the bare 
reactor and as transmitted through slab shields. 


The purpose of the ORNL-SNAP shielding program 
is the development of the technology for the nuclear 
design of radiation shields for reactor power sup- 
ply systems on spacecraft. The emphasis in the 
program is on the development and verification of 
computer programs for the design of optimum shields 
and to predict throughout the spacecraft system the 
absolute spectra of radiation resulting from reac- 
tor operation. 

Figure 1 shows a schematic diagram of the 
SNAP shielding program. The three major areas 
which are shown are cross-section development, 
transport code development and integral experiments. 
Five years ago the calculation of deep-penetration 
neutron transport in lithium hydride in a one- 
dimensional geometry was a challenging problem. 

The capabilities since that time have increased 
rapidly. 
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FIGURE l.-A Schematic Description of the ORNL-SNAP 
Shielding Program. 


The significant steps in this development in- 
clude (1) the ability to calculate neutron and sec- 
ondary gamma-ray transport in a coupled mode, (2) 
the development of the DOT (refs. 1 and 2) discrete 
ordinates code for deep-penetration two-dimensional 
transport calculations, (3) the development of a 
shield optimization program, ASOP (ref. 3), util- 
izing one-dimensional transport calculations, and 
(4) the development of the multigroup Monte Carlo 
MORSE (ref. 4) and the coupling of MORSE and two- 
dimensional DOT calculations for the analysis of 
difficult two- and three-dimensional problems. 

The current development work includes MORSE- 
DOT coupling and the optimum shaping of shields 
which is a rudimentary form of multidimensional op- 
timization. The verification of the shield design 
and analysis methods is effected through comparison 
with integral experiments. Two types of integral 
experiments have been performed at the ORNL Tower 
Shielding Facility. In the first type of experi- 
ment the measurements are of absolute gamma-ray 
pulse-height spectra produced from thin slabs which 
are exposed to a filtered neutron beam. In compari- 
son with this experiment, the calculation of trans- 
port effects is minimized and the emphasis is 
placed on the calculation of the production of sec- 
ondary gamma rays (refs. 5 and 6). 

The second type of exeriment is more prototyp- 
ical in that absolute neutron and gamma-ray spectra 
transmitted from slab shields placed adjacent to 
the TSF-SNAP reactor are compared. These compari- 
sons test the transport calculation of neutrons and 
gamma rays , as well as the production of secondary 
gamma rays. In the comparison with integral exper- 
iments, it was soon verified that the available 
cross-section data for production of secondary 
gamma rays were grossly inadequate for most mater- 
ials. This situation has now been greatly improved 
in that data for production of secondary gamma rays 
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from neutron capture as a function of neutron energy 
have been produced by use of nuclear model theory 
codes (refs. 7 and 8) and by measurements for tung- 
sten (ref. 9) and uranium (ref. 10). Current work 
is directed toward the verification and improvement 
of data for secondary gamma-ray production due to 
neutron inelastic scattering. 

As the techniques developed in the basic tech- 
nology program have become established, they have 
been incorporated in an optimum-shield-design pro- 
cedure which is exercised and demonstrated through 
detailed shield design studies. Several studies 
have been performed in which the reactor type and 
power and the radiation dose constraints are speci- 
fied and an optimized shield is designed (ref. 11). 
Through these efforts, the techniques for shield 
design are being rapidly improved. 

Also of interest are parametric optimization 
studies in which the effect of reactor size, type 
and power and the radiation dose constraints on the 
optimum shield configuration and weight have been 
studied. These studies, of course, incorporate 

considerable simplification of the geometry as 
compared to the more detailed design studies for a 

specific configuration. 

Other papers by the ORNL group at this meet- 
ing will cover, in more detail, the integral exper- 
iments for secondary gamma-ray data testing (ref . 

12) , the use of coupled MORSE Monte Carlo and DOT 
calculations for three-dimensional problems (refs. 

13 and 14) , and the design of optimized shaped 
asymmetric 4ir shields (ref. 15). In this paper, 
the status of the calculational comparisons with 
the TSF-SNAP integral experiments will be summar- 
ized. 

The source for these experiments is the TSF- 
SNAP reactor shown before assembly in figure 2. 

The reactor is a modified SNAP-2 which was designed 
and constructed by Atomics International. The 
pressure vessel shown on the left is 9 in. in diam- 
eter, 16 in. long, and contains the uranium-zircon- 
ium-hydride fuel in 36 rods with a central stainless 
steel rod containing a small quaiitity of boron car- 
bide. The reactor is reflected radially by beryl- 
lium and has four reflector control vanes. The 
coolant is NaK and the heat is removed by natural 
circulation to a NaK-to-air heat exchanger above 
the reactor. In the dry critical runs during the 



FIGURE 2. -The TSF-SNAP Reactor Prior to Assembly. 

start of the reactor, the internal power distribu- 
tion was measured by gamma-ray scans of the fuel 
elements ,by copper, uranium, and gold foil activa- 
tion and by activation of copper strips. Compari- 
sons of the power distribution deduced from these 
measurements with 05R-Monte Carlo calculations and 
DOT discrete ordinates calculations show very good 
agreement for the radial and axial power distribu- 
tions (ref. 16). 

After the assembly of the reactor was completed 
and in its final operating stage, the first measure- 
ments which were made were of the neutron angular 
current leaking from the bottom of the assembly. 
Figure 3 shows the SNAP reactor without the heat 
shield suspended from its boom. The large collima- 
tor tank contains a special 2-in. by 2-in. square 
collimator which was designed to view only 10% of 
the area of the core. By moving the collimator 
along the track and varying the angle position, 
several measurements were made of the angle-depen- 
dent neutron leakage. Data from this experiment 
showed that except for those measurements where the 
collimator was vertical and viewing the outer extrem- 
ity of the core and reflector, the spectral shape 
and intensity which were measured were very similar 
for each run. Comparisons of the measured absolute 
spectra with very detailed 05R calculations showed 
substantially good agreement (ref. 17). 
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FIGURE 3.-TSF-SNAP Reactor With NE-213 Detector 
Collimator Used For Measurement of the Angle- 
Dependent Neutron Leakage. 
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FIGURE 4 .-Comparison of Calculated and Measured 

Radial Neutron Dose Profiles Below the SNAP Reactor. 


Figure 4 shows comparison of measurements and 
calculations for lateral neutron dose traverses 
below the SNAP reactor. The measurements were made 
with a small Hornyak button, and the calculations 
were made with the two-dimensional DOT code for the 
reactor coupled, in one case with the SPACETRAN ray- 
tracing code to obtain the dose at the specified 
points, and in the other case the coupling at a 
plane just below the reactor was made to the MORSE 
Monte Carlo code which then determined the dose at 
the specified points. The agreement between the 
calculations and experiment is surprisingly good. 

In the DOT-SPACETRAN coupling, however, some ray 
effects can be observed for large radii on the plane 
close to the bottom of the reactor. 

Additional neutron and gamma-ray spectral mea- 
surements for the bare reactor are included in the 
slab shield transmission series. The neutron and 
gamma-ray slab transmission measurements have been 
performed with large slab shields supported such 
that they are located directly under the SNAP reac- 
tor. Measurements have been made with the NE-213 
fast-neutron scintillator and the sodium-iodide 
gamma-ray detector. 


Figure 5 shows a list of the experiments for 
which calculations have been performed. There are 
many other experimental runs which have been made 
but not yet calculated. From figure 5 you can see 
that neutron measurements and calculations have been 
performed for the bare reactor, and lithium-hydride, 
depleted uranium, lead, tungsten powder, heavimet 
alloy, and laminated lithium -hydride lead-uranium 
slabs. More difficult gamma-ray measurements and 
calculations have been performed for the bare reac- 
tor, and lead, depleted uranium, heavimet alloy, 
tungsten powder, and laminated lithium-hydride and 
uranium slabs. 

Figure 6 shows a scale drawing of the geometric 
configuration for the fast-neutron spectral measure- 
ments. The NE-213 detector is located deep in a 
room beneath the concrete pad . The area of the slab 
which it views is defined by the inner wall of the 
water tank. In the calculations, the leakage angu- 
lar flux from a DOT calculation of the reactor is 
used as a boundary source for the calculation of the 

slab shield. The leakage angular flux from the 
bottom of the slab shield is then integrated with 
SPACETRAN code to obtain the absolute neutron flux 
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at the NE-213 detector. In the experiment, the 
pulse-height distribution from the NE-213 is 
unfolded with the FERDOR code and compared 
with the calculation at the detector point. The 
comparisons are always on an absolute basis with 
the reactor power monitored by foils during each 
run. 

EXPERIMENTS WHICH HAVE BEEN CALCULATED 

NEUTRON 

BARE REACTOR 
6" LiH 
12" LiH 
1.5" U-238 a 
4.5" U-238 
1.5" Pb 
3" Pb 
6"Wb 
12" W 

2" HEAVIMET* 

1.5" U-238 - 6" LiH 
6" LiH -1.5" U-238 
6" LiH -1.5" Pb-1. 

6" LiH -4. 5" U-238 
6" LiH -1.5" Pb 
3" Pb-6" LiH 

a U-238 is depleted uranium, 
tyv is tungsten powder {p~ 6 gnt /cm. \ 

C HEAVIMET is tungsten-copper-nickel alloy. 

FIGURE 5. -A List of TSF-SNAP Slab Transmission 
Experiments Which Have Been Calculated. 


EXPERIMENT CONFIGURATION FOR NEUTRON MEASUREMENTS 
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DETECTOR 


FIGURE 6. -Experimental Configuration for Fast- 
Neutron Measurements. 


Figure 7 shows a comparison of experiment with 
DOT-SPACETRAN calculation and Monte Carlo 05R cal- 
culation for the fast-neutron spectra for the bare 
core. The agreement is very good except that the 
calculations are somewhat high below 1-1/2 MeV. 

Figure 8 shows a comparison of DOT calcula- 
tion with a recent measurement of the neutron leak- 
age from the bare core using the Benjamin spectro- 
meter which uses spherical hydrogen-filled detec- 
tors. The measurements were made at approximately 
36 in. below the reactor core. The calculation 
uses DOT coupled with SPACETRAN in the usual manner, 
and the agreement in intensity is quite good except 
that the calculation is somewhat high between 600 
keV and 1-1/2 MeV. Also shown is the result of 
scaling the NE-213 measurements to the hydrogen 
counter location. The agreement between the NE-213 
spectrometer and the Benjamin spectrometer in the 
overlap energy range is excellent. 



FIGURE 7 . -Comparison of Calculated and Measured 
Fast-Neutron Spectra From the Bare TSF-SNAP 
Reactor . 
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FIGURE 8. -Comparison of Calculated and Measured 
Intermediate Neutron Spectra From the Bare TSF- 
SNAP Reactor. 


GAMMA 

BARE REACTOR 
2" Pb 
3" Pb 
6" Pb 
3.5" U-238 
4. 5" U-238 
2" HEAVIMET 
6" W 

6" LiH -1.5" U-238 


5" U-238 
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Figure 9 shows a comparison of measurement 
with DOT and 05R calculations for a 12-in. -thick 
slab of lithium hydride. The overall agreement is 
very good; however, the calculations are somewhat 
high below 2 MeV and somewhat low at the higher 
energies. The multigroup lithium cross sections 
used in the DOT calculation were processed prior to 
the availability of lithium in the ENDF/B library 
but used the same data as used for the ENDF evalua- 
tion. It is important that all of the neutron- 
producing reactions in lithium be treated in order 
to obtain reasonably good agreement. In this mul- 
tigroup set, the secondary neutron distributions 
from the nonelastic reactions are assumed to be 
isotropic in the laboratory system. However, in 
the Monte Carlo calculation the angular distribu- 
tions of the secondary neutrons are treated explic- 
itly. Previous comparisons in one-dimensional 
sphere geometry between 05R and multigroup ANISN 
calculations have shown excellent agreement in 
lithium hydride to a depth of 90 cm (ref. 18). 



FIGURE 9. -Comparison of Measured and Calculated 

Fast-Neutron Spectra From a 12-in. Lithium Hydride 
Slab Under the SNAP Reactor. 


Figure 10 shows a comparison of experiment and 
DOT calculation for the fast-neutron spectra trans- 
mitted through a 2-in. slab of heavimet alloy. The 
agreement for this case is the best for any of the 
metal slabs which have been calculated. For lead, 
the calculated spectra are consistently 15% low, 
and for uranium at large slab thicknesses the calcu- 
lations are a factor of 2 to 2-1/2 low. This result 
for uranium is anomolous since comparisons of DOT 
calculations with experiment using the same slabs in 
a beam source geometry have previously shown good 
agreement when the fission multiplication in the 
slab is treated (ref. 1, page 144). Experiments are 
currently in progress to attempt to determine the 
reason for this disagreement in this case. 

Figure 11 shows a scale drawing of the geometry 
for the gamma-ray measurements and calculations. 

The 5-in. sodium iodide detector is also located 
deep in the room beneath the concrete pad . In the 
gamma-ray measurements a 10-in. thickness of borated 
polyethylene is in the beam at all times to reduce 
the neutron background in the detector . 
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FIGURE 10. -Comparison of Calculated and Measured 
Fast-Neutron Spectra From a 2-in. Slab of Heavimet 
Under the SNAP Reactor. 

EXPERIMENT CONFIGURATION FOR GAMMA RAY MEASUREMENTS 



FIGURE 11 .-Experimental Configuration for Gamma -Ray 
Measurements . 
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Calculations have been performed with DOT using 
a 27-neutron group - 60-gamma -ray group coupled 
library. Some of the calculations use the angular 
flux distribution beneath the reactor as a fixed 
boundary source for the slab shield calculations. 

In other cases, it has been necessary to calculate 
the reactor and shield as a single unit and include 
the effect of the shield on the power distribution 
in the reactor. The angular flux from the lower 
surface of the shield is integrated with the SPACE- 
TRAN code to determine the flux at the detector in- 
corporating the uncollided flux attenuation of the 
borated polyethylene. The gamma-ray flux incident 
on the detector is folded with the detector response 
function so that the calculated pulse-height distri- 
bution is compared directly with the measurement. 

Figure 12 shows the comparison of calculated and 
measured absolute pulse-height distribution for the 
bare core. The agreement below 4-1/2 MeV is excel- 
lent. Above 4-1/2 MeV, the experiment is high by 
almost a factor of 2. This disagreement is at pre- 
sent believed to be due to epithermal capture in 
structural materials in the reactor which are not 
adequately described. A comparison of dose or total 
energy flux would agree very well since the inten- 
sity is much greater at the lower energies. 



FIGURE 12. -Comparison of Calculated and Measured 
Gamma-Ray Pulse-Height Spectra From the Bare SNAP 
Reactor. 


Figure 13 shows a comparison for experiment 
and calculation of the absolute pulse-height spectra 
for a 2— in. slab of heavimet placed directly under 
the SNAP reactor. Over most of the energy range 
the calculation is 20—25% low. At the higher ener- 
gies above 7 MeV, the calculation is a factor of 2 
low. However, the integral energy flux of the cal- 
culation and experiment agrees to within 15% due to 
the fact that the calculation agrees best with the 
experiment below 2 MeV where the intensity is the 
highest. This is the best comparison for any of 
the metal slabs. 



FIGURE 13. -Comparison of Calculated and Measured 
Gamma-Ray Pulse-Height Spectra From a 2-in. 
Heavimet Slab Under the SNAP Reactor. 

The secondary gamma-ray production matrix used 
in this calculation, which included energy-dependent 
yields or epithermal capture, has been checked with 
the TSR-2 secondary gamma-ray experiment and shown 
to give substantially good results (ref. 19). 

For the other gamma-ray shields the situation 
is not nearly so good. For lead, the TSR-2 gamma- 
ray experiments have shown that the secondary 
gamma-ray production data for neutron inelastic 
scattering is entirely inadequate. For uranium, 
the comparisons of calculation and measurements for 
the SNAP reactor experiments show that the calcula- 
tions are low by about a factor of 2. This disa- 
greement may be related to the anomaly in the neutron 
spectral comparisons since the TSR-2 gamma-ray ex- 
periment shows good agreement for 238 U foils. 
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In conclusion, the comparisons and calcula- 
tions with the TSF-SNAP experiments have shown that 
the transport codes can accurately predict the neu- 
tron and gamma-ray radiation from the bare reactor. 
Also, the prediction of neutron transmission through 
shield slabs of lithium hydride, lead, and tungsten 
are adequately calculated. For secondary gamma-ray 
comparisons, the present status is not satisfactory 
in that relatively good agreement is only obtained 
for the heavimet alloy and tungsten powder assem- 
blies. 


The work presently in progress may ^oon re- 
solve the difficulties experienced in neutron and 
gamma-ray comparisons for depleted uranium slabs, 
and the disagreement for lead constitutes a low 
priority problem in that lead has been ruled out as 
a likely candidate material in an actual SNAP shield. 

Future experiments using the TSF-SNAP reactor 
will be for the verification of shield optimization 
by comparing measured and calculated dose deriva- 
tives and for verification of shield shaping pro- 
cedures. 
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